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Abstract. We consider the problem of computing the Euler characteristic of an abstract sim- 
phcial complex given by its vertices and facets. We show that this problem is #P-complete and 
present two new practical algorithms for computing Euler characteristic. The two new algorithms 
are derived using combinatorial commutative algebra and we also give a second description of them 
that requires no algebra. We present experiments showing that the two new algorithms can be 
implemented to be faster than previous Euler characteristic implementations by a large margin. 



1. Introduction 

The Euler characteristic of a topological space is an invariant used in a variety of contexts 
such as category theory, algebraic geometry and differential geometry. In combinatorics, the Euler 
characteristic of a simplicial complex is related to the Mobius function of a poset and the inclusion- 
exclusion principle |16] and to valuations on simplicial complexes [12 to name but a few connections. 

The reduced Euler characteristic of an abstract simplicial complcsQ A is 

X (A) - J2 = -/-I + ./o - ./i + /2 - /3 + • • • 

o-eA 

where fi denotes the number of faces (elements) of dimension i in the complex^ The dimension of 

a face a is dim(cr) |cr| — 1. 

In Section [2] we prove that computing the Euler characteristic of a simplicial complex specified 
by its vertices and facets is #P-complete, which is a formal way of stating that Euler characteristic 
is a difficult computational problem. We also show that the problem of deciding if x (A) = is 
not in NP unless #P is no harder than NP. This answers two open questions posed by Kaibel and 
Pfetsch in their survey [TT] . 

In Section[3]we introduce two new practical algorithms for computing Euler characteristic. These 
two algorithms were conceived of in terms of combinatorial commutative algebra, and Section [3] is 
written solely in terms of algebra. Section |3] independently describes the same two algorithms in 
terms of simplicial complexes and without any reference to algebra. Section [S] describes how the al- 
gebra was translated to simplicial complexes and how doing so brought up interesting mathematics. 



Date: December 21, 2011. 

^An abstract simplicial complex A is a family of sets closed under taking subset, so if tr G A and r C c then 
r € A. This is closely related to the notion of a simplicial complex which is a set of polyhedra with certain properties. 
All complexes in this paper are abstract. See Section |4. II for further background. 

^The usual definition of Euler characteristic is x(^) /o — /l /2 — /s -I- ■ ■ ■ • The difference is that x{^) does 
not count the empty set, while X (^) does, so X (A) = x('^) ~ 1- A.11 Euler characteristics in this paper are x(^) 
rather than x(^) because that simplifies the formulas. 
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Finally, Section [5] presents experiments that show that the two new algorithms can be implemented 
to be faster than previous Euler characteristic implementations by a large margin. 

2. The Complexity of Euler Characteristic 

We describe the complexity class #P and then prove that Euler characteristic is #P-complete. 
This is a precise way of saying that Euler characteristic is a difficult computational problem. We 
also consider the complexity of decision problems associated to Euler characteristic. See Section 
14.11 for basic definitions relating to simplicial complexes. 

The complexity of Euler characteristic has been studied before, but not of a simplicial complex 
specified by its vertices and facets. It has been studied for the case of the input being a CW-complex 
specified as a circuit [6] in the context of real valued computation and for the input being a sheaf 
[T] in the context of algebraic geometry. 

2.1. The complexity class #P. The complexity class #P is the set of counting problems asso- 
ciated to decision problems in NP. For example the decision problem "does a logical formula have 
some satisfying assignment of truth-values?" is in NP, while "/low many satisfying assignments of 
truth-values does a logical formula have?" is in #P. The former is called SAT while the latter is 
called #SAT. A problem is #P-complete if it is in #P and any other problem in #P can be reduced 
to it in polynomial time. 

There is already a list of problems that are known to be #P-complete, which is very helpful 
when proving that a new problem is #P-complete, as then a problem in ^P is #P-complete if some 
other #P-complete problem reduces to it. For example it is known that #SAT is #P-complete even 
when restricted to formulas with two literals per clause and no negations [18\ A SAT formula is a 
conjunction of clauses, where each clause is a disjunction of some number of literals. For example 

(aV-.6) A(aVc)A(-.6Vc), 

where a, b and c are boolean variables. Here the satisfying truth assignments (a, b, c) are 

{(0,0,1), (1,0,1),(1,1,0), (1,1,1)}. 

The output for SAT with this formula as input is "yes" since there is a satisfying truth assignment. 
The output for ^Skl is "4" , since there are four satisfying truth assignments. The output for #SAT 
does not include the satisfying truth assignments themselves, only the number of them. 

The program for the rest of this section is to formally define a problem EulerChar in #P that 
represents the Euler characteristic problem, and then to prove that EulerChar is #P-complete. 

2.2. Euler Characteristic is in #P. The most straigtforward way to define EulerChar would 
be to have the input be the facets and vertices of a simplicial complex A and have the output be 
simply X (^)- It is immediate that this could never be in because x (^) can be negative while 
#P is a class of counting problems so that their output must be a natural number. 

To arrive at a satisfactory definition of EulerChar, the first step is to observe that 

(1) X (A) XI (-l)'*""^"^ = # (odd faces) - # (even faces) . 

Read # (odd faces) as "the number of odd faces of A" , where a set is odd if it has an odd number 
of elements which is to say that its dimension is even. It is not hard to argue that counting the 
number of even faces is a problem in #P, and that counting the number of odd faces is a problem 
in =f/=P as well. Unfortunately, we know of no theorem stating that a difference of two functions in 
#P is again in #P. So we must find an alternative way to express the Euler characteristic. 
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Let A have n vertices. Then 

# (even faces) + # (even non- faces) = # (even sets) = 2"^"^, 

which together with Equation (jlj imphes that 

X (A) + 2"""^ = # (odd faces) + # (even non-faces) . 

Consider the decision problem "does A have an odd face or an even non-face?" . This problem 
is in NP where a certificate of a "yes" -answer is any concrete odd face or even non- face. Define 
EulerChar to be the counting version of this. The input is then the vertices and facets of A and 
the output is the number of odd faces and even non-faces, that is the output is x (^) + 2"^^. We 
can subtract 2"~^ in polynomial time, which justifies that EulerChar represents the problem of 
computing the Euler characteristic of a simplicial complex. 

We conclude that computing Euler characteristic is a problem in #P when expressed formally in 
the form of the EulerChar problem. 

2.3. Euler Characteristic is #P-complete. The main result in this section is Theorem [T] which 
states that EulerChar is ^P-complete. This is an example of the fact that even trivial problems 
can have a counting version that is #P-complete. To see that the problem "does A have an odd 
face or an even non-face" is especially trivial, observe that only A = fails to have the even face 0. 

Theorem 1. EulerChar is #P-complete. That is, the problem of computing the Euler characte- 
ristic of a simplicial complex given by its vertices and facets is #P-complete. 

Proof. We proved in Section [2?2] that EulerChar is in #P. We prove the statement of the theorem by 
showing that the ^P-coniplete problem ^SAT reduces to EulerChar. We introduce an intermediate 
problem IndepSum and prove that #SAT reduces to IndepSum and then that IndepSum reduces to 
EulerChar. Given a SAT formula S, the combination of these two reductions yields a simplicial 
complex A such that x (A) is the number of truth assignments that satisfy S. 

We need to introduce some terminology. Let the parity sum of a set of sets S be P{S) 
^ggg(— 1)1*1. For example the parity sum of a simplicial complex A is P{A) — — x (A). Let G be 
a simple graph with vertex set V. Then a set of vertices S* C is dependent if it contains both 
endpoints of some edge of G. Otherwise S is independent. 

We can now define the problem IndepSum. The input of IndepSum is the vertices and edges of a 
graph G, and the output is the parity sum of the set of independent sets of G. 

IndepSum reduces to EulerChar: Let G be a simple graph with vertex set V and define a 
simplicial complex A such that the facets of A are the complements of the edges of G. Then a set 
of vertices is a face of A if and only if the complement contains an edge, that is if and only if the 
complement is a dependent set of G. 

If a C_ V then let a V \ a be its complement. Let D be the set of dependent sets of G or 
equivalently I? = {ct|(t e A}. Since (-1)'^' = (-l)"(-l)l'^l and ~x{A) is the parity sum of the 
faces of A, we get that 

P(i?)-(-l)"P(A) = -(-l)"x(A). 

Let / be the set of independent sets. Every set is either dependent or independent, so assuming 
that y 7^ we get that {V {V) is the set of all subsets of V) 

P{I) + P{D) ^ P{I UD)=^ P{r {V)) = 0. 

We conclude that P{I) -P{D) = (-l)"x (A) so that IndepSum reduces to EulerChar. 
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#SAT reduces to IndepSum: Let S he a SAT formula. We construct a graph G such that the 
number of truth assignments that satisfy S equals (—1)" times the parity sum of the independent 
sets of G. 

Let vi, . . . ,Vn be the variables that appear in the formula S and let ci, . . . ,Ck be the clauses 
that appear in S. For each variable Vi we introduce a 3-clique with vertices Ti, Fi and Di. Here Tj 
represents Vi having the value true and Fi represents false. For each clause Cj we introduce a vertex 
Cj . If the literal Vi appears in clause Cj with no negation, then we add an edge between Ti and Cj . 
If the literal -^Vi appears in clause Cj then we add an edge between Fi and Cj. We claim that the 
number of truth assignments that satisfy S equals (—1)" times the parity sum of the independent 
sets of this graph G. 

For concreteness, consider the SAT formula 

{vi V -1U2) A {vi V Vz) A (-1U2 V U3). 

The graph that we construct based on this formula is shown in Figure [TJ 

Let A be the set of vertices named Di or Cj and let B be the set of vertices named Ti or Fj . Let 
/ be the set of independent sets of C and let Ib be the set of independent sets that are subsets of 
B. Define the function p: I ^ Ib hy p{d) d\A. 

We are going to prove that i) if p"^(c?) = {d} then \d\ = n and ii) that the set of such d is in 
bijection with the set of truth assignments that satisfy S. We are also going to prove Hi) that if 
p^^{d) 7^ {d} then the parity sum P{p^^{d)) is zero. These three statements imply that 

# (satisfying truth assignments) = (-1)" ^ P{p-^{d)) = (-1)"P(/), 

dels 

where we use that is a partition of /. It only remains to prove i), ii) and Hi). 

i) If p~^{d) — {d} then \d\ — n: Suppose that d G Ib such that p~^{d) — {d}. Pick some 
variable Vi. Then d U {Di} is dependent since otherwise it would be an element of p~^{d). As Di 
is only adjacent to T,; and Fi, it must be the case that d contains one of Ti and Fi. It cannot 
contain both as there is an edge between them. If d contains Ti then we assign the value true to 
Vi and otherwise d contains Fi and we assign the value false to Vi. In this way d encodes a truth 
assignment to the variables of the formula S. 

ii) {d E lb \p~^id) — {d}} is in bijection with the satisfying truth assignments of S: 
Pick some clause Cj. Then dU {Cj} is dependent so d must contain some Ti or Fi that is adjacent 
to Cj and this implies that the truth assignment that d represents satisfies the clause Cj. This 




Ci:Viy^V2 C2-ViVV3 C3-^V2\/V3 

Figure 1 . Illustration for the proof of Theorem [1] 
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establishes a bijection between the set oi d E Ib such that p ^{d) = {d} and the set of truth 
assignments that satisfy S. 

in) ii p-^{d) ^ {d} then P(p-^(d)) — 0: Let d e h such that p~^{d) =1 {d}. Then we can 
pick some vertex a G A such that dL) {a} £ p~'^{d). Then d does not contain any vertex that is 
adjacent to a, and there are no edges between the elements of A, so if we let 

E = {h e p-\d) \a(^h}, f"^' {he p-\d) \aeh] 

then h ^ h\J {a} is a bijection from i? to F so that P{F) = —P{E). As {E, F} is a partition of 
p-^id) we then get that P{p-'^{d)) = P{E) + P{F) = 0. □ 

The Euler characteristic is the alternating sum of the entries of the /-vector, so Euler characte- 
ristic reduces to /-vector. So we get the following result of Kaibel and Pfetsch |Tl] as a corollary. 

Corollary 2. The problem of computing the /-vector of a simplicial complex given by its vertices 
and facets is #P-hard. 

2.4. Decision Problems. In this section we investigate the complexity of decision problems as- 
sociated to Euler characteristic. Kaibel and Pfetsch pose the open problem of whether deciding 
X (A) = lies in NP [11]. Theorem |3] answers this question in the negative unless #P is no harder 
than NP. It is a central conjecture of computational complexity theory that #P is harder than NP. 

Theorem 3. Let Eq be the problem of deciding if x (A) = where A is a simplicial complex given 
by its facets and vertices. Then Eq is co-NP-hard. Also, Eq does not lie in NP unless #P is no harder 
than NP. 

Proof. Let E< be as in Lemma [5] and assume that Eq is in NP. Then E< is in NP n co-NP by 
Lemma [5] This allows us to compute Euler characteristic in NP n co-NP using binary search. Euler 
characteristic is =i^P-complete by Theorem [l] so then =i^P is no harder than NP n co-NP. 

Eq is co-NP-hard: Let 5 be a SAT formula. The proof of Theorem [1] constructs a simplicial 
complex A such that x (^) is the number of satisfying truth assignments to S. So the NP-complete 
problem SAT reduces to the decision problem x (^) 7^ 0. So the complement of Eq is NP-hard, which 
implies that Eq is co-NP-hard. □ 

This leaves an open problem of whether Eq lies in co-NP, since Theorem [3] does not rule that 
out. If Eq does lie in co-NP, it would then be proven that NP ^ co-NP unless ^'P is no harder than 
NP since Eq would then lie in co-NP and not in NP. It is an open problem whether NP 7^ co-NP. 

Theorem 4. Let E>o be the problem of deciding if x (^) > where A is a simplicial complex 
given by its facets and vertices. Then E>o is #P-hard. 

Proof. Let E> and E< be as in LemmaO The argument used to prove the equivalence of Eo and E= 
in Lemma [5] also works to show that E>o and E> are equivalent. Then in particular both E> and 
E< reduce to E>Oj so Euler characteristic reduces to E>o using binary search. Euler characteristic 
is #P-complete by Theorem [T] so then E>o is #P-hard. □ 

Lemma 5. Consider the following decision problems, where A is a simplicial complex given by its 
facets and vertices and k is an integer, 

Eo: x(A) = 0, E<: x (A) < /c, 

E=: x(A) = fc, E>: x{^) > k. 

If any one of these problems are in NP then they are all in NP n co-NP. There are polynomial time 
reductions in both directions between Eq and E= and between E< and E>. 
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Proof. Eo and E= are equivalent: Assume without loss of generality that A 7^ 0. Use Lemma IH] 
to construct a polynomial size simplicial complex F such that x (r) = fc — 1 and A n F = {0}. Then 

^' A U F has X (^') = X (A) + X (F) + 1 = X (A) + fc so x (A) = if and only if x (*) = k. This 
gives a polynomial time reduction in both directions between Eq and E= and also shows that Eq is 
in NP or co-NP if and only if E= is in NP or co-NP respectively. 

E< and E> are equivalent: Let F be a simplicial complex such that x (r) — —1 and such that 
the set of vertices of F is disjoint from the set of vertices of A. Let =' A © F as in Theorem [TSl 
whereby x (^) = —X {^)- Then x (A) < fc if and only if x (^) > —k. This gives a polynomial time 
reduction in both directions between E< and E> and also shows that E< is in NP or co-NP if and 
only if E> is in NP or co-NP respectively. 

E= in NP E< in NP: If E= is in NP then we can certify the exact value of x (A), which will 
also serve as a certificate for x (A) > k. If E< is in NP, then so is E> in which case we can certify 
that fc — 1 < X (A) < k + 1 which serves as a certificate of x (A) = k. 

If one problem is in NP, then they are all in co-NP: Assume that one of Eq, E=, E< and 
E> is in NP. Then we have shown that they are all in NP. So we know that E= is in NP, which 
allows us to certify the exact value of x (A). This also serves as a certificate for when x (A) < fc is 
not true and when x (A) = fc is not true, so E< and E= are in co-NP. We have already proven that 
this implies that E> and Eq are also in co-NP. □ 

Lemma [6] constructs a simplicial complex with a given Euler characteristic k such that the bit 
size of the complex is bounded by a fixed polynomial in the bit size of the Euler characteristic which 
is [log fc] . It is necessary to bound the bit size of the simplicial complex in this way since otherwise 
the proof of Lemma [5] would not go through. 

Lemma 6. Let k be an integer. Then there is a simplicial complex A such that x (A) — k and A 
has no more facets and no more vertices than 2P + 31 + 7 where I — [log2(|fc|)] or Z = if fc = 0. 

Proof. The proof is based on inclusion-exclusion along with Theorem 1151 
The case A; = 0, 1: Let Aq = and Ai ({^i} , {^2}). 

The case k — 2": Let n be a positive integer and define F^ {{xnii} , {xni2} , {xnis}) ■ Let 
= Fi © ■ • • ® F„. We have x {^^) = 2 so X (*n) = 2". Also has 3n vertices and 3n facets. 

The case fc > 0: Write k — (hi- - -60)2 in binary such that k — X]L=o^"2"' ^" ^ ^^"^ 
hi^Q. This implies that / ~ [log2(fc)]. 

Let W he a. finite set of non-empty simplicial complexes with disjoint vertex sets. Then the 
only face in more than one element of W is 0, so x (UVK) = X]yiew^(^) + 1^1 ^ 1- 



Then x = (p — 1) * ( — 1) = ~ by Theorem [TSl Observe that $ has no more facets and no 
more vertices than \W'^\ + 4, < I + 4,. Let Wn U {$} and A^ UW„. Then 




{*„ \hn = 1} we have x {^^W^) = k+ jW/J - 1. Let p ^ |VF/J + 1 and 

{{y^} ...,{yp})®{W,b},{a,c},{b,c}). 





n=0 
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The case A; < 0: Let fl {{zi, Z2} , {zi, z^} , {z2, z^}) and observe that x (^) = ^1- Let 
Afc © A_fc so that x(A/j) = — x(A_fc) = k. Observe that A^ has no more vertices and no 
more facets than 

l + ^{1(1 + l)) + 3=^f + ^l + 7 < 2f + 31 + 7. □ 

3. Algebraic Algorithms for Euler Characterisic 

In this section we describe two new algorithms for computing Euler characteristic of a simplicial 
complex using algebraic techniques. In Section |4] we present these same two algorithms in the 
language of simplicial complexes. This section is independent from Section |4] and it only uses 
algebra. 

More precisely the two algorithms we present in this section compute the coefhcient of x = 
xi ■ ■ ■ Xn in the multigraded Hilbert-Poincare series numerator H (/) of a square free monomial 
ideal /. In Section [5] we show that this is equivalent to computing the Euler characteristic of a 
simplicial complex. For that reason we define 

X (I) coefficient of x in H (/). 

The summary of what Section [5] shows in detail is that given a simplicial complex A we can define 
a monomial ideal / such that X (/) = X (^)- Computing / from A takes little time. See Section [S] 
for details on the relationship between the algebraic algorithms in this section and the simplicial 
algorithms in Sectional 

3.1. Background and Notation. We work in a polynomial ring k[xi, . . . , a;„] over a field k and 
with variables xi, . . . , x„. A monomial ideal is a polynomial ideal generated by monomials. Let 
/ be a monomial ideal. Each monomial ideal has a unique minimal set of monomial generators 
min (/). The exponent vector of a monomial m is a vector v such that m — YVi=i ■ ^ monomial 
has full support if it is divisible by x. A monomial ideal has full support if lcm(min(/)) has full 
support. The colon of two monomials a and 6 is a : 6 icm(^a,fe) ^ rpj^^ colon of a monomial ideal / 
by a monomial a is 

I : a — {m \ma G / } = (m : a |to G / and m is a monomial) . 
The multigraded Hilbert-Poincare series hilb (/) is the possibly infinite sum of standard mono- 
mials of /, that is hilb (/) X^m^/ where the sum is taken over monic monomials m. The 
multigraded Hilbert-Poincare series can be written as a rational function 

hilb (/) = 

^ ' {l-Xi)---{l~Xn) 

where H (/) is a polynomial called the multigraded Hilbert-Poincare series numerator. 

3.2. Divide... Both algorithms we present are divide-and-conquer algorithms. They take a mono- 
mial ideal / and split it into two simpler monomial ideals J and K such that X (I) ^ X (J) + X (K)- 
This process proceeds recursively until all the remaining ideals are simple enough that they can be 
processed directly. 

Let J? be a square free monomial and let / be a square free monomial ideal. The divide steps for 
the two algorithms are derived from the equation 

(2) hilb (/) = (hilb {I ■.p))p + hilb (/ +{p)). 
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By giving these three Hilbert-Poincare series the same denominator (1 — a;i) •••(!— a;„), we get a 
similar equation for the Hilbert-Poincare series numerators 

H(/) = (H(/:p))p + H(/+(p)). 

By considering the coefficient of x on both sides, we then get that 

x{i)^x{{i-p)p) + x{i + {p))- 

It would simplify this expression if we could write x{I : p) instead of x ((/ : p)p). This does not 
work directly since in general x{I : p) will be zero since / : p will not have full support. We are 
working within a polynomial ring R = ] . We will argue that if we change the ring that 

I : p is embedded in, then it becomes true that xi^ ■ P) = xii^ '■ p)p)- 

Since the variables that divide p do not appear in min (/ : p) we can embed I' I : p into a 
ring R' =^ k[P] where P =' {xi \xi does not divide p}. Let x' HP = ^ be the product of the 
variables in R' . Then x{I') — x{{I ■ p)p) so we consider / : p to be embedded in i?', which gives 
us the final equation behind splitting 

(3) X{I) = X{I-P)+X{I+^))- 

There are two different algorithms for Hilbert-Poincare series that are based on Equation We 
present two analogous algorithms for Euler characteristic that are based on Equation ([3]). 

The Hilbert-Poincare series algorithm due to Bigatti, Conti, Robbiano and Traverso [SJ |3] uses 
Equation ([2|) directly as we have written it. It is a divide-and-conquer algorithm that splits a 
monomial ideal / into the two simpler monomial ideals / : p and I + {p). We call this algorithm the 
BCRT algorithm for Hilbert-Poincare series. In analogy with that algorithm, we propose a BCRT 
algorithm for Euler characteristic that uses Equation ([3]) directly as written - it splits / into the two 
simpler ideals / : p and / -I- (p). We call it the algebraic BCRT algorithm for Euler characteristic. 

We call the p in Equation ([3]) the pivot. Section l375l explores strategies for selecting pivots. 

There is also a Hilbert-Poincare series algorithm due to Dave Bayer and Michael Stillman [3] 
that we will call the DBMS algorithm for Hilbert-Poincare series. It is based on writing Equation 
© as 

hilb (/ +{p))= hilb (/) - hilb {I -.p). 

Given an ideal J, the idea is to choose p e min (J) and let / =' (min (J) \ {p}) such that 

hilb (J) = hilb (/ +{p))^ hilb (/) - hilb {I :p). 

In this way J splits into the two simpler ideals / and I . p. In the same way, we can rewrite 
Equation ^ as 

(4) xiJ)^x{I+{p))^x{I)~xiI-P)- 

We propose a DBMS algorithm for Euler characteristic that uses this equation to split J into / and 
/ : p. We call it the algebraic DBMS algorithm for Euler characteristic. Note that the pivots in the 
DBMS algorithm are minimal generators of the ideal, which would not make sense for the BCRT 
algorithm. 

3.3. ... and Conquer. A square free monomial ideal / is a base case for both algorithms when 
Theorem [7] applies. Note that the improvements in Section [3761 enables further base cases. 

Theorem 7. Let J be a square free monomial ideal. Then 
(1) if / does not have full support then x (/) = 0, 
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(2) if / has full support and the minimal generators min (/) of / are pairwise prime monomials 
then x(/) = 

(3) if / has full support and |min (/)| —2 then x {I) = 1- 

Proof. (HI): All the monomials with non-zero coefficient in H (/) can be written as lcm(M) for some 
M C min (/). If / does not have full support then neither does any monomial of the form lcm(Af). 
Then x must have a zero coefficient since it has full support. 

(12): If the elements of min (/) are relatively prime then H (/) = nm6min(/) (-'^ ~ "^)- ideal 
has full support we then get that x = n„gmi„(/)m so the coefficient of x is (— . 

([2]): If {a, b} =' min (/) and g gcd(a, b) then I ^ g (c, d) where c'=^ a : g and d'=^ b : g. As c 
and d are relatively prime by construction, we get that 

H (/) = <?H ((c, d)) = g{l - c)(l - d) = cdg -eg- dg + g. 

So the coefficient of x = cdg in H (/) is 1. □ 

These base cases improve on the ones for the Hilbert-Poincare series algorithms in that they 
apply more often and can be processed more quickly. For example the Hilbert-Poincare series 
algorithms as well have a base case when the elements of min (/) are relatively prime, but it takes 
exponential time to process that base case since nmGmin(7) (-'^ ~ '^^^^ have 2l'°'"^^^l terms. Here 
all that is required is to determine if |min (/)| is even or odd. The base case when / does not have 
full support does not exist for the Hilbert-Poincare series algorithms. 

3.4. Termination and Complexity. It is clear that the DBMS algorithm terminates since |min (/) 
decreases strictly at each step. For the BCRT algorithm, termination requires that we choose the 
pivots p such that 1 ^ p ^ I since otherwise we get an infinite number of steps from I — I : p or 
I = I+{p)- 

If we cannot choose a p such that 1 ^ p ^ I then I = (1) which is a base case. If 1 ^ p ^ / and / 
has full support, then I C. I + (p) and I C I : p. So if the BCRT algorithm does not terminate then 
there would be an infinite sequence of strictly increasing ideals in contradiction to the fact that the 
ambient polynomial ring is Noetherian. So both algorithms terminate. 

We have seen that the DBMS algorithm gets rid of at least one minimal generator at each step. 
The BCRT algorithm gets rid of at least one variable at each step if the pivot is chosen to be a 
single variable Xi. It is immediate that Xi is not a variable of / : x^. To see that Xi can also be 
removed from I + (xi), observe that the only minimal generator that is divisible by Xi is Xi itself, so 
Xi is (/ -l- (xi))-independent from the other variables and so can be removed using the independent 
variables technique from Section [3.61 

Since the Euler characteristic problem is ^P-complete we expect both algorithms to run in at least 
single exponential time. Using the transpose technique from Section 13.61 we can interchange the 
number of variables n with the number of minimal generators |min (/)|. So if I == min(n, |min (/)|) 
then both the BCRT and DBMS algorithms have 0{q2^) asymptotic worst case time complexity 
where g is a polynomial. So both algorithms run in single exponential time. We expect that a more 
careful analysis could reduce the base of the exponential. 

3.5. Pivot Selection. We have proven that the BCRT and DBMS algorithms terminate in a finite 
amount of time. To be useful in practice, the amount of time until termination should be small 
rather than just finite. The strategy used for selecting pivots when splitting an ideal has a significant 
impact on performance. We describe several different pivot selection strategies here and compare 
them empirically in Section [51 
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A popular variable is a variable that divides a maximum number of minimal generators of the 
ideal. In other words, a popular variable Xi maximizes |min (/) n {xi)\. A rare variable is a variable 
Xi that minimizes |min (J) n {xi)\ with the constraint that Xi ^ min (/). 

If there are several candidate pivots that fit a given pivot selection strategy, then the pivot used 
is chosen in an arbitrary deterministic way among the tied candidates. 

BCRT Pivot Selection Strategies. Recall that BCRT pivots p satisfy 1 ^ p ^ I. 

popvar: The pivot is a popular variable, 
rarevar: The pivot is a rare variable. 

random: The pivot is a random variable e such that e ^ min (/). 

popgcd: Let xi be a popular variable. The pivot is the gcd of three minimal generators 
chosen uniformly at random among those minimal generators that Xi divides. 

The strategies popvar and popgcd have been found to work well for the BCRT algorithm for 
Hilbert-Poincare series, so we also try them here, rarevar and random have been included to have 
something to compare to. 

DBMS Pivot Selection Strategies. Recall that DBMS pivots are elements of min(/). 

rarevar: The pivot is a minimal generator divisible by a rare variable. 

popvar: The pivot is a minimal generator divisible by a popular variable. 

msLXSupp: The pivot is a minimal generator with maximum support. 

minsupp: The pivot is a minimal generator with minimum support. 

random: The pivot is a minimal generator chosen uniformly at random. 

rarest: The pivot is a generator that is divisible by a maximum number of rare variables. 
Break ties by picking the generator that is divisible by the maximum number of second- 
most-rare variables and so on. 

raremax: The pivot is chosen according to rarevar where ties are broken according to 
maxsupp. 

3.6. Improvements. We present several improvements to the DBMS and BCRT algorithms. 

Independent Variables. We say that two subsets A, B C_ {xi, . . . , x„} are I -independent if min (/) 
is the disjoint union of min (/ H k[A]) and min (/ H ^[i?]). This is another way of saying that the 
minimal generators of / can be partitioned into two subsets such that every generator in one set is 
relatively prime to every generator from the other set. If A and B are /-independent then 

hiib (/) = hiib (/ n k[A]) ■ hiib (/ n k[b]) . 

This is a standard technique for computing Hilbert-Poincare series [3l O I15j . and it applies to 
computing Euler characteristic as well since 

x{i) = xiin4A])-x{inK[B]). 

The existence of an /-independent pair {A, B) can be determined in nearly linear time |15) . and 
such a pair can cut down on computation time dramatically. However, independence rarely occurs 
for large random ideals and detecting it does take some time, so this technique is not worth it unless 
there is some reason to suspect that it will apply to a given ideal. 



COMPLEXITY AND ALGORITHMS FOR EULER CHARACTERISTIC OF SIMPLICIAL COMPLEXES 11 



Eliminate Unique Variables. Suppose Xi divides only one minimal generator p G min(J). Use 
Equation (|4]) with p as the pivot to get that 

x(J) = x(/+(p» = x(/)-x(/:p) = -x(/:p), 
since x (!) = as no generator of / is divisible by Xi so / does not have full support. 

Transpose Ideals. Let M be a matrix whose entries are or 1. Each row of M is then a 0-1 vector 
that we can interpret as the exponent vector of a square free monomial. Define (M) to be the 
monomial ideal generated by the monomials whose exponent vectors are the rows of AI. On the 
other hand, given a square free monomial ideal /, we can take the exponent vectors of the minimal 
generators of / and put them in a matrix. Define Mj to be that matrix. We label the rows of 
Mi with the elements of min (/) and we label the columns of Mj with the variables in the ambient 
polynomial ring of /. Let the transpose ideal /"^ of / be the ideal generated by the transpose of the 
matrix of /, so trans (/) == (M/"^). 

Theorem |8] states that x{I) ^ X (trans (/)), so we can transpose the ideal without changing the 
Euler characteristic. The BCRT algorithm is more sensitive to the number of variables than it is 
to the number of generators, so it can be beneficial to transpose the ideal if it has fewer generators 
than variables. The DBMS algorithm is opposite of this in that it is more sensitive to the number 
of generators than the number of variables. 

Another situation where transposing can be beneficial is in the case where a column of Mj 
dominates another column. If we take the transpose, those two columns will become generators 
and the dominating generator will not be minimal. When we then take the transpose again, we 
will have fewer variables than we started with. This process can repeat itself several times as the 
removal of dominating columns from the matrix can cause rows to be dominating, and removing 
those dominating rows can then cause yet more columns to become dominating. 

Theorem 8. If / is a square free monomial ideal then x{I) = X (trans (/)). 

Proof. The proof is by induction on the number of variables n. Let / be a square free monomial 
ideal. Choose a variable Xi and let 

J (min (/) \ (xi)) = (m G min (/) \xi does not divide m) . 
The plan of the proof is to show that 

(5) X{I)^X{I■.X^)-X{J■■X^) 

and that 

(6) X (trans (/)) = x (trans (/ : Xi)) — x (trans (J : Xi)) . 

Recall that we embed / : Xi and J : Xi in a subring that does not have the variable Xi, so the result 
follows from these two equations by applying the induction assumption to / : Xi and J : xi. It only 
remains to prove Equation ([5]) and Equation ([5]). 

Equation ([5]): Equation [3] with the pivot implies that 

XiI)^X{I ■.Xi)+xiI+{x^)). 

Now J does not have full support and I + (xi) = J + (xi) , so we get by Equation (jl]) that 

X{I+ (a;,)) = X ( J + {x,)) = X ( J) - X ( J : x,) ^ -x{J : x,) . 

Equation This part of the proof is easier to follow by the reader drawing pictures of the 
matrices involved. Let v be column Xi of Mj. Then the entries of v are indexed by min (/) and if 
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m G min (/) then Vm = 1 if and only if Xi\m. Let A be the result of removing column Xi from Mj. 
Then (M/^) = (^4^) + (x") so Equation g]) implies that 

X (trans (/)) = x ((m/)) = x ((A^) + (xl) = X ((A^)) - X ((^^) : ^'') ■ 

The colon / : Xi corresponds to removing column Xi of Mj so (A) — I : Xi. The colon can also 
reduce the number of minimal generators, so Mj.^xi can have fewer rows than A does. However, 
those extra rows are exponent vectors of non- minimal generators so they do not impact {A) . Then 
Lemma ini implies that 

X (trans (/ : Xi)) = x (trans ((A))) = x {{A^)) ■ 

It now suffices to prove that trans {J : Xi) = {A^) : x". Let B be the result of removing column 
Xi of Ml and also those rows m 6 min (/) such that = 1. We see that Mj is M/ with the same 
rows removed as for B. Doing a colon by Xi removes column so J : = (i?). 

It remains to prove that (i?^) = {A^) : x". Observe that {A^^ : x"" removes those columns 
m £ min (/) of where = 1 which are the same columns that are removed from B"^ . Both 
and B^ have had row Xi removed so {^) '■ x" = {B^). We have proven that 

trans (J : x^) = trans ((B)) = {B^) = {A^) : x" 
where trans {{B)) = (S^^ depends on the observation that no row of B dominates any other. □ 
Lemma 9. If (A) = {B) then x ((^^)) = X ((-S^)) i^i matrices A and B. 

Proof. Assume without loss of generality that no row of B dominates any other. Then A has all 
the rows that B does and also some additional non-minimal rows. Assume by induction that there 
is only one additional row r. Let d be some other row of A that is dominated by r. Then d and r 
contribute variables Vd and Vr to the ambient ring of (A^). We get by Equation [3] that 

x{{A^))^x{{A'^)--Vr)+x{{A^) + {vr)). 

All generators of (^"^) that are divisible by Vd are also divisible by Vr, so (A^) + (wr) does not 
have full support at Vr so x {{A'^) + (vr)) = 0. Therefore we have that 

X{{A-))^X{{A^)--V^)^X{{B^)), 
where we are using that the colon {A'^'j : Vr corresponds to removing row r from A. □ 

Base Case for |min (/)| = 3. Assume that all unique variables have been eliminated, that / has full 
support and that |min(/)| ~ 3. Then every variable Xi divides 2 or 3 elements of min(/). We can 
ignore the variables that divide all 3 minimal generators as they make no difference to the Euler 
characteristic. For every minimal generator there must be at least one variable that does not divide 
it. So after removing repeated variables by taking the transpose twice, we see that / must have the 
same Euler characteristic as (xy, xz, yz). So x (I) ^ 2. 

Partial Base Case for |min(/)| — 4. Suppose that |min(/)| — 4, that every variable divides exactly 
two elements of min (/) and that the number of variables is 4. Then x (^) — ~1- There should be 
more rules like this, though identifying them by hand is laborious and error prone. 

Make a Table. It would be beneficial for each small k to perform a computer search to make a 
table of all ideals / with |min (/)| = k up to reordering of the variables and the various techniques 
for simplifying an ideal that we have presented. Then the Euler characteristic of ideals with few 
generators could be computed through a table look-up. 
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Data Structures. The exponents of / are all or 1, so we can pack 32 or 64 exponents into a single 
32 or 64 bit machine word, and in that representation many operations become much faster. We 
used this standard technique in our implementation and while the general concept is simple we 
warn that the implementation details are tricky. 

For sparse or complement-of-sparse exponent vectors, it might pay off to only record the zero 
entries and one entries respectively, though this is not something that we have pursued. 

4. SiMPLiciAL Algorithms for Euler Characteristic 

In this section we present two algorithms for Euler characteristic that work directly with simplicial 
complexes. These two algorithms are equivalent to the monomial ideal based algorithms from 
Section [3l We introduce the simplicial versions from the ground up, so this section is independent 
from Section [3] and does not use any algebra. See Section [S] for more on the connection between 
the algebraic and simplicial versions of the two algorithms. 

4.1. Background and Notation. In the introduction we wrote that a simplicial complex is a 
finite set of sets that is closed under subset. We are going to set up an algebra-simplicial dictionary, 
and for that to work we associate a set of vertices to a simplicial complex. This way a simplicial 
complex can have vertices that are not an element of any of its faces. Section 15.11 shows why this 
is necessary. Definition 1101 adds this vertex set structure to simplicial complexes. 

Let V {V) be the set of all subsets of V for any set V. 

Definition 10. Given a finite set V, a simplicial complex A is a subset of P {V) that is closed 
under taking subsets. In other words, A is a set of subsets of V such that if cr G A and t C a then 
T S A. The elements of A are called faces and the elements of Va V are called vertices of A. 

The facets of A are the maximal faces of A with respect to inclusion. The set of facets is 
denoted by fac(A). The closure (D) of a set of sets D C V (V) is the smallest simplicial complex 

that contains F, namely (F) UdeoP (d). A simplicial complex is uniquely given by its facets 
since A = (fac (A)). The complement of a set a C Va is Va \ cr. Two sets a, r C Va are 
co-disjoint if their complements are disjoint, or equivalently if a U r = Va. If t C Va then 

A e T = {cr e A |a n r = 0} = {cr \ r |a e A} . 

The vertex set of A r is Vagt =^ Va \ t- 

We will use that Euler characteristic respects inclusion-exclusion in the sense that 

X (A U A') = X (A) + X (A') - X (A n A') . 

4.2. Divide... Both algorithms we present are divide-and-conquer algorithms. They take a simpli- 
cial complex A and split it into two simpler complexes A' and A" such that x (A) — x {^')+X (A")- 
This process proceeds recursively until all the remaining complexes are simple enough that they 
can be processed directly. Splitting (the divide step) is based on Theorem [TTl 

Theorem 11. If is a non-empty set of vertices then 

X(A) =x(Aea)+x(AUP(a)). 

Proof. First observe that 

AnV (ct) = {t e A\t C a} = {t e A\Tna ^9} ^ AQa. 
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□ 



In analogy with the algebraic BCRT algorithm from Section [3l we can use the equation in 
Theorem [TT] as written to split a simplicial complex A into the simpler complexes A ct and 
A yP (a). We will refer to this algorithm as the simplicial BCRT algorithm. It is not immediately 
clear that A QW and AuV (a) are simpler than A is. For now, consider that A a has fewer 
vertices if cr C Va and that AUV (cr) has fewer non- faces if ct ^ A. 

We call the a in Theorem [Til the pivot. Section l475l explores strategies for selecting pivots. 

We can also write the equation in Theorem [TT] as 



In analogy with the algebraic DBMS algorithm from Section [31 we can use this equation to split a 
simplicial complex D into the simpler complexes A and A ct. We will refer to this algorithm as 
the simplicial DBMS algorithm. Note that the pivots in the simplicial DBMS algorithm are facets 
of the complex, which would not make sense for the simplicial BCRT algorithm. 

4.3. ... and Conquer. A simplicial complex A is a base case for both algorithms when Theorem 
[T2| applies. The improvements in Section [476] enable further base cases. 

Theorem 12. Let A 7^ be a simplicial complex. Then 

(1) if A is a cone then x (A) = 0, 

(2) if A is not a cone and the facets of A are pairwise co-disjoint then x (A) = (— 

(3) if A is not a cone and |fac (/)| = 2 then x (A) — 1. 

Proof. ([T]): This is well known and not hard to prove. 
By induction on |fac(A)| using Lemma [T3l 
([3|): The two facets are co-disjoint when ignoring unused vertices. □ 

Lemma 13. Let D he a, simplicial complex and let a S fac {D) such that a ^ Vd is co-disjoint to 
every other facet of D. Then x (D) = — x (A a) where A =^ (fac (D) \ {cr}). 

Proof. Let v ^ a and r G fac (A). Then u G r so A is a cone and x (A) — 0. By Theorem [TT] 



4.4. Termination and Complexity. It is clear that the simplicial DBMS algorithm terminates 
since |fac(A)| decreases strictly at each step. For the simplicial BCRT algorithm, termination 
requires that we choose the pivots a such that a C Va and cr ^ A since otherwise we get an infinite 
number of steps from A = A0aorA = AuP (cr) 

If we cannot choose a cr such that cr C Va and cr ^ A then A = P (Va) which is a base case. If 
we always choose pivots a such that a C Va and ct ^ A then termination of the simplicial BCRT 
algorithm follows from the fact that either the number of vertices |Va| or the number of non- faces 




X(i?) =x(AU7'(p)) = x(A)-x(A0ct) 



x(A0ct). 
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|A| — 2l^'^l — |A| decreases between any two steps. In fact some consideration shows that |A| 
decreases strictly at each step. 

We have seen that the DBMS algorithm gets rid of at least one facet at each step. The BCRT 
algorithm gets rid of at least one vertex at each step if the pivot is chosen to be of the form a == {e}. 
It is immediate that e is not a vertex of A Bo". To see that e can also be removed from AUV (cr), 
observe that a is the only facet that does not contain e, so e is (A U P ((T))-independent from the 
other vertices and so can be removed using the independent vertices technique from Section 14.61 

Since the Euler characteristic problem is #P-complete we expect both algorithms to run in at 
least single exponential time. Taking nerves as described in Section 14. 6i we can interchange the 
number of vertices |Va| with the number of facets |fac (A)|. So if I min(|VA| , |fac (A)|) then both 
the BCRT and DBMS algorithms have 0{q2') asymptotic worst case time complexity where q is 
a polynomial. So both algorithms run in single exponential time. We expect that a more careful 
analysis could reduce the base of the exponential. 

4.5. Pivot Selection. We have proven that the BCRT and DBMS algorithms terminate in a finite 
amount of time. To be useful in practice, the amount of time until termination should be small 
rather than just finite. The strategy used for selecting pivots when splitting an ideal has a significant 
impact on performance. We describe several different pivot selection strategies here and compare 
them empirically in Section |6l 

A rare vertex is a vertex that belongs to a minimum number of facets. A popular vertex is a 
vertex e that is an element of a maximum number of facets with the constraint that {e} ^ fac(A). 

If there are several candidate pivots that fit a given pivot selection strategy, then the pivot used 
is chosen in an arbitrary deterministic way among the tied candidates. 

BCRT Pivot Selection Strategies. Recall that BCRT pivots a satisfy a C Va and <t ^ A. 

The reason that some of these names seem opposite of their definition is that the names come 
from the algebraic setting and the translation to simplicial complexes involves taking a complement. 
For example a rare variable of an ideal translates to a popular vertex of the corresponding simplicial 
complex. 

popvar: The pivot is {e} for e a rare vertex, 
rarevar: The pivot is {e} for e a popular vertex. 

random: The pivot is {e} for e a random vertex such that {e} ^ fac (A), 
popgcd: Let e be a rare vertex. The pivot is the union of three facets chosen uniformly at 
random among those facets that do not contain e. 

The simplicial BCRT algorithm presented here is based on the algebraic BCRT algorithm for 
Hilbert-Poincare series, and the strategies popvar and popgcd have been found to work well for 
Hilbert-Poincare series computation. So we also try them here, rarevar and random have been 
included to have something to compare to. 

DBMS Pivot Selection Strategies. Recall that DBMS pivots are elements of fac (A) . 

rarevar: The pivot is a facet that does not contain some popular vertex, 
popvar: The pivot is a facet that does not contain some rare vertex, 
maxsupp: The pivot is a facet of minimum size, 
minsupp: The pivot is a facet of maximum size, 
random: The pivot is a facet chosen uniformly at random. 
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rarest: The pivot is a facet that lacks (does not contain) a maximum number of popular 
vertices. Ties are broken by picking the facet that lacks the maximum number of second- 
most-popular vertices and so on. 

raremax: The pivot is chosen according to rarevar where ties are broken according to 
maxsupp. 

4.6. Improvements. We present several improvements to the DBMS and BCRT algorithms. 

Independent Vertices. If A and F arc simplicial complexes with disjoint vertex sets Va and Vr, then 

A © r = (A X P (Vr)) U {P (Va) x T). 

So if A = {{x} , {y}) and T = {{a} , {b} , {c}) then 

A ® r = {{x, a, b, c} , {y, a, b, c} , {x, y, a} , {x, y, b} , {x, y, c}) . 

Using Theorem [T5] we can compute % (A © F) in terms of x (^) ^nd x (r). So if we are computing 
X (^) for a simplicial complex ^, then we could significantly simplify the task by finding simplicial 
complexes A and F such that = A © F. 

We say that two subsets AjBCV^ are 'i> -independent if fac (^P) is the disjoint union of 

Fa {cr e fac {^)\a\A^A} and Fb = {cr e fac {^i) \a\B = B} . 

If \I/ = A © F then it is immediate that Va and Vr are ^'-independent. On the other hand, if A and 
B are \l/-independent, then = A © F where 

A = (Fa) © a and F (Fb) © B. 
We have proven Theorem I14[ which together with Theorem 1151 generalizes Lemma [T3l 

Theorem 14. A simplicial complex can be written as ^ = A © F if and only if there is a 
^'-independent pair {A,B). 

Theorem 15. Let A and F be simplicial complexes with disjoint non-empty vertex sets Va and 
Vr. Then 

x(A©F) = x(A)x(r). 

Proof. Choose any set D such that A = {D) . Then we get by inclusion-exclusion that 

X (A) = X ( u ^ (^) ) - E ( n ^ i'^yj 

\cr£D I vCD \<Jev / 

vCD 

Let Fa = fac (A) x {Vr} and Ft = {Va} x fac (F). Then A © F = (i^A U Fr) so in the same way 
we get by inclusion-exclusion that 

x(A©F)= {~^P^'x{v{nv)) 

= E E(-i)'"^'''^'x(^(n(dug))). 

dCF^ gCFr 

Let d' C fac (A) and g' C fac (F) such that d^d'UVr and g ^ g' UVa- We will prove that 

(8) X (r (n(d u g))) = -X (P {nd')) x (P (ng')) • 
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Since |d| — and \g\ — \g'\ we then get that 

x(A©r)= E (-i)''''+''''x(^(nd'))x(7'(n5')) 

d'Cfac(A) 3'Cfac(r) 

- E (-i)l''l+'x(^(nrf')) E (-i)''''^'x(7'(n5')) 

d'Cfac(A) g'Cfac(r) 

= x(A)x(r). 

It only remains to prove Equation ([8]). Observe that x (''')) = if t is non-empty and otherwise 
X (t)) = So Equation ([5]) is equivalent to the statement that 

n((i U .g) = ^ nrf' and n 5' = 0. 

As every face in d contains Vr and every face in g contains Va, the only way for n((i U to be 
empty is if nd — Vr and = Va- The statement then follows from the observation that flc? = Vy 
if and only if C\d! — and likewise ng = Va if and only if fig' = 0. □ 

Eliminate Abundant Vertices. Suppose e is an element of every facet of D except a G fac (D). Use 
Equation ([7]) with a as the pivot to get that 

Xp) =x(AUT'(a)) =x(A)-x(Aea) =x(Aea), 

since x (A) = as every facet of A contains e so A is a cone. 

Take Nerves. The nerve of a simplicial complex A is 

nerve (A) = {v C fac (A) |nw / } . 

A complex and its nerve have the same homotopy type [ini Thm. 10], so x (A) — x (nerve (A)). 
Taking the nerve corresponds to transposing the facet-vertex incidence matrix M. Taking the nerve 
twice will remove dominated columns from M and then remove any dominated rows that might 
have appeared. This process can continue as removing rows of M can cause yet more columns to 
be dominated. 

If no columns are dominated then taking the nerve once will exchange the number of vertices and 
facets. The BCRT algorithm is more sensitive to the number of vertices than it is to the number 
of facets, so it can be beneficial to do computations on the nerve if the complex has fewer facets 
than vertices. The DBMS algorithm is opposite of this in that it is more sensitive to the number 
of facets than the number of vertices. 

Base Case for |fac (A)| — 3. Assume that all abundant vertices have been removed, that A is not 
a cone and that |fac (A)| ~ 3. Then every vertex that actually occurs in the complex is an element 
of one facet. The nerve of the nerve of A will then have the form ({a} , {b} , {c}). So x (A) = 2. 

Partial Base Case for |fac (A)| — 4. Suppose that |fac (A)| — 4, that every vertex is an element of 
exactly two facets and that the number of vertices is 4. Then x (A) = — 1. There should be more 
rules like this, though identifying them by hand is laborious and error prone. 

Make a Table. It would be beneficial for each small k to perform a computer search to make a 
table of all simplicial complexes A with |fac(A)| — k up to reordering of the vertices and the 
various techniques for simplifying a simplicial complex that we have presented. Then the Euler 
characteristic of simplicial complexes with few facets could be computed through a table look-up. 
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Data Structures. In our implementation we have represented a facet as a 0-1 vector, where we pack 
32 or 64 entries into a single 32 or 64 bit word. Many operations are very fast in this representation. 
While the general concept is simple we warn that the implementation details are tricky. 

For sparse or complement-of-sparse vectors, it might pay off to only record the zero entries or 
one entries respectively, though this is not something that we have pursued. 



5. Ideals to Simplicial Complexes 

In this section we describe in more detail how the simplicial algorithms in Section|4]are equivalent 
to the algebraic algorithms in Section [3l We treat this topic in detail in part to support our claim 
that the algebraic and simplicial algorithms are equivalent and in part because the algebra-simplicial 
translation brings up interesting mathematics. 

To set up the algebra-simplicial translation, let A be a simplicial complex with vertices xi, . . . ,Xn 
that are simultaneously the variables in a polynomial ring. The algebra-simplicial connection is via 
the function (f> from sets to monomials defined by (f> (cr) =^ IIct. In other words, (j) (a) is the product 
of variables (vertices) that are not in a. So if n = 3 then cf) (0) = 2:1X2X3 and (j) ({xi, X3}) = X2. 

Extend </) from sets to simplicial complexes by 0(A) =^ ('/'(c) |cr € A) . Then is a bijection 
from the facets of A to the minimal generators of (p (A) , so (j) (A) can be quickly computed given 
A. We can also describe (f>{A) as the Stanley- Reisner ideal of the Alexander dual of A, though we 
will not use this alternative description of (A) here. 

The fundamental observation that we have been using is that x (4' (^)) — X (^) so that we can 
compute X (^) with monomial ideal methods on (j) (A). 

Theorem 16. If A is a simplicial complex then x (^) = X (</> i^)), where x {4> (^)) is the coefficient 
of x in the multivariate Hilbert-Poincare series numerator H (/) of /. 

Proof. The proof is based on a formula of Bayer concerning upper Koszul simplicial complexes. 

Let m be a monomial and let / be a monomial ideal. Then the upper Koszul simplicial compZeall 
[m [2] of / at TO is 

' x™ 

(T C {Xl, . . . ,X„} ^ e ^ 

The numerator of the multivariate Hilbert-Poincare series of / is related to the Euler characteristics 
of the Koszul complexes of / by the following formula due to Bayer [2j Proposition 3.2] 



, def 



H(/)-l= ^ X(A^. 



Since A^*'^^ — A, we get that 



Coefficient of x in H (</. (A)) = x [^t^^^) = X (A) . □ 

5.1. The Inverse of 0. In Section HTT] we modified the definition of a simplicial complex to have 
an associated set of vertices so that a vertex need not actually appear in any of the faces of the 
complex. Here we give an example that shows that this change is necessary for to be a bijection 
between simplicial complexes and square free monomial ideals. In short, the vertex set definition 
is necessary to preserve information about the variables in the ambient polynomial ring. 



'This notion of Koszul complex is not to be confused with the chain complex notion of a Koszul complex. 
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We have two functions (jj, one that maps sets to square free monomials and one that maps 
ideals to square free monomial ideals. They are defined by (j) (v) '= Ilv and (f> (A) (0 {v) |f e A) 
respectively. The inverse of A iH- (/) (A) is given by 

0-i(/) = {„C{a;i,...,x„}|0(w)G/}. 

We define the vertex set of cj)'^ (/) to be the set of variables in the ambient ring of /, even if 
some of those variables do not appear in any face of (/). We give an example that shows that 
A I— (A) would not have an inverse if the vertex set of (j>~^ (/) were the union of its faces. 

Consider the ideal / =' (a;ia;2, xi^a) C k[xi,X2,X3] and let A (/) so that fac (A) = 

{{xs} , {x2}}- The question now is what the vertex set of A is. The union of faces is {x2, xs} so if 
that were the set of vertices then (jj ({2:2}) = X3 and ({3:3}) — X2 so then (j) (/)) = (j) (A) = 
{x3,X2) ^ /. We observe that if the vertex set were the union of faces, then there would be no 
simplicial complex A such that (j) (^) = I- So (f> would be a bijection not from simplicial complexes 
to square free monomial ideals, but from simplicial complexes to square free monomial ideals with 
the property that gcd(min (J)) — 1. 

Using the definition that we have given, the vertex set of A is still {xi,X2,X3} even though xi 
does not appear in any face of A. With our definition (j) does have as an inverse, and indeed 
we see that ^ ((^"^ (/))=</> (A) = /. 

Even if the input ideal / to the algebraic BCRT and DBMS algorithms has the property that 
gcd(min(/)) = 1, the intermediate ideals that these algorithms generate do not necessarily have 
that property. So if we had defined the vertex set to be the union of faces, then we would have had 
to deal with this issue in some other way. 

5.2. Divide... The algebra algorithms are based on Equation [3] which states that for monomials p 

xii) = xii ■P) + Xil+{P))- 
The simplicial algorithms are based on Theorem [11] which states that for sets a 

X(A) =x(Aea)+x(AUP(a)). 
These two equations are equivalent since if A = (/) and cr = (p) then 

5.3. ... and Conquer. Theorem [7] specifies base cases for the algebra algorithms when 

(1) / does not have full support, 

(2) / has full support and the minimal generators min (/) are pairwise prime, 

(3) / has full support and |min(/)| = 2. 

Theorem [T^] specifies base cases for the simplicial algorithms when 

(1) A is a cone, 

(2) A is not a cone and the facets of A are pairwise co-disjoint, 

(3) A is not a cone and |fac (/)| = 2. 

If A ~ (j)^^ (/) then the algebra and simplicial conditions are equivalent. To be more precise 

(1) / has full support if and only if A is not a cone, 

(2) monomials a and b are relatively prime if and only if <p^^ (a) and cf)^^ (b) are co-disjoint, 

(3) |min(/)| = |fac(A)|. 
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5.4. Termination and Complexity. The algebraic DBMS algorithm terminates as |min(7)| de- 
creases strictly at each step. The simplicial DBMS algorithm terminates as |fac(A)| decreases 
strictly at each step. These reasons are equivalent when A = (/) as then |min(/)| = |fac(A)|. 

The algebraic BCRT algorithm terminates if the pivots p are chosen such that p ^ 1 and p ^ I. 

The simplicial BCRT algorithm terminates if the pivots a are chosen such that a C Va and cr ^ A. 
If A = (f>^^ (I) and a = (p) then these conditions are equivalent. 

5.5. Pivot Selection. The pivot selection strategies that have the same name in the two sections 
on pivot selection are equivalent. The main points in proving this are that for a monomial ideal / 
and a simplicial complex A =^ (7) 

• a variable Xi is rare if and only if the vertex Xi is popular, 

• a variable Xi is popular if and only if the vertex Xi is rare, 

• a minimal generator m has maximum support if and only if (/>~^ (m) has minimum size, 

• a minimal generator m has minimum support if and only if (m) has maximum size. 

5.6. Improvements. Let / be a square free monomial ideal and let A =^ (^~^ (/). The two sections 
on improvements present equivalent techniques in the same order. 

Independent Variables — Independent Vertices. Let A and B be disjoint sets of variables. Let 
I C k[A] and J C k[B]. Then I + J C k[A U iJ] and the algebraic section on independent 
variables proves that x(7 + J) = x{I)x{J)- Furthermore, an ideal K can be written as a sum 
7 + J if and only if there arc disjoint sets A and B such that min (K) is the union of min (K) r]K[A] 
and min (K) n k[B]. In that case we say that A and B are 7C-independent and then TsT = 7 + J for 
7 = (min (K) n k[A]) and J = (min (K) n 

For the simplicial side of things, let A and F be simplicial complexes such that the vertex sets 
Va and Vr are disjoint. Then 

A © r = (A X (Vr)) U {V (Va) x F) and ^Aer = Va U Vr 

The simplicial section on independent vertices proves that x (A © F) = x (A) x (F). Furthermore, 
an ideal ^ can be written as A © F if and only if there are disjoint sets A and B such that fac {^) 
is the disjoint union of 

Fa = {a e fac (*) \a\A = A} and Fb = {a € fac (*) \a\B = B}. 

In that case * = A © F for A = {Fa) and F = (Fb). 

Addition of ideals and © of simplicial complexes are equivalent. If A (j}"^ (7) and F =^ (/)~^ (J) 
then A © F = (7 + J). Observe that we are using three different functions here as the 
ambient polynomial ring is different for each of them. Recall that the definition of involves 
taking a complement, and if u is a set of variables then the meaning of the complement v depends 
on what the variables in the ambient polynomial ring are. 

Independence of variables and vertices are also equivalent. If =^ (K) and A and B are sets 
of variables/vertices, then A and B arc Ti'-indepcndent if and only if A and B are ^'-independent. 

Let A and B be TC-indepcndent and let 5' = (f>^^ [K] so that A and B are also ^I'-indepcndent. 

Let A = {Fa) and F = {Fb). Let 7 = (min (7C) n and J = (min (7C) n Then 

Ti" = 7 + J and * = A © F. Furthermore, A = ^"^ (7) and F = (J). 

The notion of 7-independence is standard, though we have not been able to find any reference 
in the literature to ^'-independence or the operation © on simplicial complexes. 
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Eliminate Unique Variables — Eliminate Abundant Vertices. A variable Xi is unique in / if and 
only if Xi is an abundant vertex of A. 

Transpose Ideals — Take Nerves. Both the nerve and the transpose of an ideal are transposing a 
matrix where the columns are variables/vertices and the rows are generators/facets. One of these 
matrices can be derived from the other by replacing all O's by 1 and simultaneously replacing all 
I's by 0. This shows that (j>~^ (trans (/)) = nerve (A). 

As far as we know, the transpose operation has not been applied to monomial ideals before - we 
are investigating if a monomial ideal and its transpose have any interesting relations between them 
when the ideal is not square free. 

Base Case for |min(/)| = 3 — Base Case for |fac(A)| = 3. These are equivalent as |min(/)| — 
|fac (A) I and both base cases give the same Euler characteristic of 2. 

Partial Base Case for |min (/)| =4 — Partial Base Case for |fac (A)| = 4. These are equivalent. 
Make a Table. This is the same idea. 

Data Structures. These are the same considerations. Sparse exponent vectors correspond to large 
facets, while small facets correspond to mostly-one exponent vectors. In particular, it is not the 
case that sparse exponent vectors correspond to small facets. 

6. EMPIRICAL Evaluation of Euler Characteristic Algorithms 

We have implemented the BCRT and DBMS algorithms for Euler characteristic in the program 
Frobby [14|. In this section we explore the pivot selection strategies for both algorithms, and we 
compare the BCRT and DBMS implementations to several other systems with Euler characteristic 
functionality. 

We use simplicial terminology in this section. To recover equivalent statements using algebraic 
terminology read "monomial ideal" for "complex" , read "variables" for "vertices" , read "minimal 
generators" for "facets", read "having not full support" for "being a cone" and read "transpose" 
for "nerve" . 

We compare the following implementations, listed in alphabetical order. 

Frobby version 0.9.0 [14j : Frobby is a free and open source system for computations on 
monomial ideals. We have written a C-I--I- implementation of the BCRT and DBMS algo- 
rithms for Euler characteristic in Frobby. We ran Frobby with the option -minimal turned 
on in Table [6] as all faces given in the input are facets. 

GAP version 4.4.12 [9]: GAP is a free and open source system for computational discrete al- 
gebra. It computes the Euler characteristic of a complex by enumerating all faces. The im- 
plementation is written in the GAP scripting language. We use the SCEulerCharacteristic 
function from version 1.4.0 of the simpcomp package [7]. We extend the memory limit for 
GAP to 2 GB. 

Macaulay 2 version 1.4 [8]: Macaulay 2 is a free and open source computer algebra system. 
It computes the Euler characteristic of a simplicial complex using the BCRT algorithm 
for Hilbert-Poincare series to get the /-vector. The time consuming parts of the code are 
written in CH — h. Due to some inefhciencies in the SimplicialComplexes package it is faster 
to compute the Euler characteristic using the poincare function directly on a monomial 
ideal instead of going through a simplicial complex. We have used this faster method. 
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Sage version 4.7 [17 : Sage is a free and open source computer algebra system. It computes 
the Euler characteristic of a simphcial complex by enumerating all faces of the complex. 
The implementation is written in Python. We use the eulerCharacteristic function from 
the SimplicialComplexes package. 

We use the following complexes for the comparison. 

random-a-6: A randomly generated complex with a vertices and h facets. The complex is 
generated one facet at a time. A prospective facet a is generated at random so that each 
vertex has a 50% chance to be an element of tr. If a is contained in or contains any previously 
generated facet then a is discarded and another prospective facet is generated in its stead. 

We generated these example using the genideal action of Frobby. 
nicgraph-a-6: The simplicial complex of all not 6-connected graphs on a given set of a vertices. 
Here we view each possible edge as a vertex of the simplicial complex. A graph is fe-connected 
if it is connected and it cannot be disconnected by removing 6—1 edges. 

We genenerated these examples with the NotlConnectedGraphs function in Sage. 
rook-a-6: The a x 6 chessboard complex. 

Let V and W be sets of vertices such that |y| = a and |T4^| = 6. A partial matching 
between V and is a set of pairs (w, w) with v & V and w £ W such that each vertex 
is in at most one pair. The set of partial matchings between V and W forms a simplicial 
complex called the chessboard complex. It corresponds to ways of placing rooks on an a x 6 
chessboard so that no rook attacks any other rook. 

We initially generated these examples with the ChessboardComplex function in Sage. 
We could not generate large enough examples with that function so we added the same 
functionality to Frobby and used that, 
match-a: The matching complex on a vertices. 

Let y be a set of vertices with \v\ = a. A partial matching on y is a set of pairs {w, w} 
for v^w such that each vertex is in at most one pair. The set of partial matchings on 
V forms a simplicial complex called the matching complex. 

We initially generated these examples with the MatchingComplex function in Sage. We 
could not generate large enough examples with that function so we added the same func- 
tionality to Frobby and used that. 

In all tables |Va| refers to the number of vertices and |fac(A)| refers to the number of facets. 
OOM stands for "out of memory" and indicates that the computation was terminated due to the 
program reporting an out of memory error. All experiments were run on an Intel(R) Core"'"^ 2 Duo 
CPU T7500 at 2.20GHz with 4GB of RAM running Mandriva Linux 2010.1. 

We checked that all the programs gave the correct answer for every input. All times are the 
median time of three runs. Time is in each case measured by the programs themselves, except for 
Frobby that was timed using the Unix time command line utility. All times exclude the time taken 
to start the program and read the input file, except for the Frobby times that do include startup 
and input time. We had to cut out startup and input time as it was taking more time than the 
Euler characteristic computation itself in some cases, and these experiments are not intended to be 
about starting programs or reading input. We did not do this for Frobby as Frobby starts up and 
reads input in a tiny fraction of the time taken for the Euler characteristic computation. 

6.1. Pivot Selection Strategies for the BCRT and DBMS Algorithms. The BCRT and 

DBMS algorithms are parameterized on a pivot selection strategy that determines which pivot to 
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I^aI 


|fac(A)| 


random 


popvar 


rarevar 


popgcd 


4000 


20 


0.05 


0.05 


0.06 


0.06 


6000 


20 


0.07 


0.07 


0.07 


0.07 


8000 


20 


0.09 


0.09 


0.08 


0.09 


10000 


20 


0.11 


0.11 


0.11 


0.12 


4000 


30 


3.70 


2.31 


5.73 


9.83 


6000 


30 


4.84 


3.25 


6.31 


13.41 


8000 


30 


5.87 


4.18 


7.61 


15.56 


10000 


30 


7.25 


5.34 


9.59 


20.53 


20 


4000 


0.13 


0.12 


0.12 


0.24 


20 


6000 


0.24 


0.24 


0.22 


0.45 


20 


8000 


0.37 


0.37 


0.36 


0.65 


20 


10000 


0.55 


0.55 


0.52 


0.95 


30 


4000 


4.11 


2.63 


5.78 


14.53 


30 


6000 


6.04 


4.06 


7.69 


22.01 


30 


8000 


7.48 


5.37 


9.11 


28.85 


30 


10000 


9.48 


6.99 


10.83 


33.41 


100 


100 


7.07 


1.38 


40.15 


166.56 


120 


120 


32.01 


5.92 


170.28 


2109.31 


140 


140 


130.44 


23.19 


740.48 


>7200 


160 


160 


491.43 


89.87 


2729.22 


>7200 


180 


180 


1497.50 


261.09 


>7200 


>7200 


200 


200 


4965.33 


796.03 


>7200 


>7200 


220 


220 


>7200 


1756.16 


>7200 


>7200 


240 


240 


>7200 


5051.55 


>7200 


>7200 



Table 1. BCRT pivot selection strategies. All times are in seconds. 



use at each step. We compare the strategies described in Section l3.5| /Section 14.51 on a battery of 
randomly generated complexes. 

BCRT Pivots. Table [T] shows that popvar is the best BCRT pivot selection strategy for these 
randomly generated complexes, popvar is already faster for the simpler ideals and its lead over the 
other strategies increases with the number of facets and especially with the number of vertices. 

random: This strategy is here to be able to tell if the other strategies are better or worse 

than a random choice of vertex, 
popvar: We believe that popvar performs well because a rare vertex e gives A U {{^}^ 

few facets as possible. 

rarevar: rarevar does the opposite of what popvar does. It is never far ahead and is mostly 
significantly behind popvar. For the balanced examples it is much worse than even a 
random choice of vertex. 

popgcd: popgcd works well for the Hilbert-Poincare series algorithms that the algorithms we 
present here are based on. However, the data shows that it is often significantly worse than 
choosing a random vertex when it comes to computing Euler characteristics. 
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DBMS Pivots. Tabic [2] shows that raremax is the best DBMS pivot selection strategy for these 
randomly generated complexes. 

random: This strategy is here to be able to tell if the other strategies are better or worse 

than a random choice of pivot, 
raremax: raremax combines the benefits of removing a facet that lacks a popular vertex with 

the benefit of removing a small facet. Table [2] shows that this combination is better than 

the pure strategies rarevar and maxsupp. 
rarest: We believe that rarest performs well because removing facets that lack popular 

vertices tends to create complexes that are close to being a cone, 
rarevar: rarevar is almost as good as rarest, which is reasonable since rarest is a more 

sophisticated way of breaking ties in rarevar. 
popvar: popvar does the opposite of what rarevar does, so it is reasonable that it does 

poorly. 

maxsupp: We believe that maxsupp performs better than random because removing small 

facets tends to create complexes that are close to being a cone, 
minsupp: minsupp performs worse than random, confirming that it is beneficial to select 

small pivots. 

6.2. Variations of the BCRT and DBMS Algorithms. Table [3] and Table H show times for 
the pivot selection strategies when the nerve technique from Section l3.6| /Section 14.61 is turned off. 
First of all we observe that turning the nerve technique off does not change the ranking of the pivot 
selection strategies. Furthermore, we see that without nerves the BCRT algorithm is much more 
sensitive to the number of vertices while the DBMS technique is much more sensitive to the number 
of facets. The nerve technique hides this sensitivity as it allows to interchange the number of facets 
and the number of vertices, so the DBMS algorithm can adjust the input to make the number of 
facets less than the number of vertices and vice versa for the BCRT algorithm. 

An alternative to the nerve technique is to use a hybrid approach where the BCRT algorithm 
is used for ideals with more facets and the DBMS algorithm is used for ideals with more vertices. 
If we compare the tables we see that the DBMS algorithm with the nerve technique turned on is 
faster than the hybrid approach even for ideals with more facets than vertices, so in this experiment 
the hybrid approach is inferior to the nerve technique. 

6.3. Comparison of Euler Characteristic Implementations. Table [6] compares several imple- 
mentations of Euler characteristic algorithms. This also serves as a comparison of the algorithms 
used by these implementations. Evaluating the practicality of an algorithm as opposed to an im- 
plementation is difficult because quality of implementation has a significant effect on performance 
yet quality of implementation cannot easily if at all be measured or corrected for. 

An example of an implementation (as opposed to algorithm) difference is that Frobby and 
Macaulay 2 are written in C-I--I- that compiles to native code while Gap and Sage do not compile 
to native code due to the languages that they are written inQ While the choice of implementation 
language can make a significant difference for performance, the magnitude of the differences in 
Table [5] is so large that choice of language is unlikely to be the main factor in our estimation. 



Sage does have modules written in Cython which is similar to Python but that does compile to native code. 
Many components of both Sage and Gap are written in native languages such as C and CH — h. However, that is not 
the case for the Euler characteristic components of Gap and Sage. Parts of Macaulay 2 are written in the interpreted 
Macaulay 2 language, but the Hilbert-Poincare series code is written in C-|--|-. 
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\Va\ 


|fac(A)| 


random 


popvar 


maxsupp 


rarest 


raremax 


minsup 


rarevar 


4000 


20 


0.07 


0.06 


0.06 


0.04 


0.04 


0.08 


0.04 


6000 


20 


0.08 


0.08 


0.07 


0.06 


0.06 


0.09 


0.05 


8000 


20 


0.09 


0.09 


0.08 


0.07 


0.07 


0.11 


0.07 


10000 


20 


0.11 


0.11 


0.10 


0.09 


0.09 


0.15 


0.10 


4000 


30 


4.53 


4.46 


2.36 


1.14 


1.06 


13.34 


1.24 


6000 


30 


5.96 


6.77 


3.58 


1.47 


1.46 


15.62 


1.63 


8000 


30 


7.83 


8.28 


4.46 


1.97 


1.90 


21.17 


2.26 


10000 


30 


10.50 


10.18 


5.64 


2.59 


2.53 


27.11 


3.24 


20 


4000 


0.19 


0.18 


0.16 


0.12 


0.12 


0.22 


0.12 


20 


6000 


0.32 


0.31 


0.29 


0.23 


0.24 


0.39 


0.24 


20 


8000 


0.48 


0.48 


0.46 


0.37 


0.37 


0.57 


0.38 


20 


10000 


().()() 


().G5 


().G3 


0.52 


0.54 


U.74 


U.34 


30 


4000 


5.08 


4.97 


3.10 


1.88 


1.76 


10.68 


2.06 


30 


6000 


7.70 


7.44 


4.81 


2.97 


2.85 


16.25 


3.29 


30 


8000 


10.16 


10.02 


6.46 


3.99 


3.98 


19.61 


4.48 


30 


10000 


12.63 


13.17 


8.90 


5.00 


4.98 


25.35 


5.58 


100 


100 


5.89 


5.84 


1.29 


1.05 


0.79 


40.83 


1.13 


120 


120 


25.42 


26.19 


5.79 


4.72 


3.44 


173.82 


5.06 


140 


140 


108.40 


119.19 


22.43 


18.43 


13.03 


749.66 


21.34 


160 


160 


427.46 


464.28 


87.13 


69.92 


49.59 


2767.93 


78.52 


180 


180 


1291.26 


1241.91 


253.15 


192.95 


143.14 


>7200 


232.63 


200 


200 


4211.65 


4137.68 


769.79 


614.39 


434.43 


>7200 


755.41 


220 


220 


>7200 


>7200 


1704.65 


1440.19 


980.31 


>7200 


1665.74 


240 


240 


>7200 


>7200 


4871.15 


4022.04 


2697.84 


>7200 


4811.50 



Table 2. DBMS pivot selection strategies. All times are in seconds. 



|Va| |fa'C(A)| random popvar rarevar popgcd 
30 10000 7^50 5^36 9^35 3024" 
10000 30 >7200 62.13 >7200 >7200 

240 240 > 7200 6 961.76 >7200 >7200 

Table 3. BCRT pivot selection strategies without nerves. All times are in seconds. 

|Va| |fac(A)| random raremax rarest rarevar popvar maxsupp minsupp 

30 10000 31811 52^39 643! 77\87 243^1 71^85 2382.40 

10000 30 10.64 2.58 2.66 3.26 11.30 5.97 25.95 

240 240 >7200 4184.07 6212.34 >7200 >7200 >7200 >7200 



Table 4. DBMS pivot selection strategies without nerves. All times are in seconds. 
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We can draw the firm conclusion from Table [5] that the BCRT and DBMS algorithms presented 
in this paper can be implemented to be faster on this set of complexes than all the other Euler 
characteristic implementations that we have compared. We have found no faster Euler characteristic 
implementations than these, so we believe that the comparison we have made is the most fair and 
informative comparison that can be made using the implementations that exist today. 

From Table m Tableland Tabled we see that the DBMS algorithm is faster than the BCRT 
algorithm for all the complexes when both algorithms use their best pivot selection strategy. So we 
can recommend using the DBMS algorithm over the BCRT algorithm, even though the difference 
is slight. The nerve technique is vital to make the DBMS algorithm always be faster - without it, 
good performance could only be reached by implementing both algorithms and choosing which to 
use based on the ratio of facets to vertices. 

The implementations in Sage and GAP are based on enumerating faces. Table[5]shows the number 
of faces of each complex, and there is a clear trend that the time taken by Sage and GAP is related 
to the number of faces of the complex. In contrast the times for the implementations in Frobby 
and Macaulay 2 do not have a clear relationship to the number of faces. 

For example nicgraph-8-2 has 2928 times more faces than match-13 does, and as expected Sage 
and GAP take much longer to compute the Euler characteristic of nicgraph-8-2 than of match-13. 
In contrast Frobby computes the Euler characteristic of nicgraph-8-2 in 0.06s while it takes 7.84s 
for match-13. So nicgraph-8-2 has 2928 times more faces than match-13 does yet it takes 260 times 
less time to compute its Euler characteristic using Frobby. We give this as evidence that the time 
taken by the BCRT and DBMS algorithms for Euler characteristic depends more on the structure 
of the complex than on the number of faces of the complex. 

We find that the performance of Frobby and Macaulay 2 in this comparison lends credence to 
the idea of using algebraic formulations and implementations for combinatorial problems. 
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Example 


Vertices 


Facets 


Faces 


X(A) 


rook-6-6 


36 


720 


13,327 


185 


rook-7-7 


49 


5,040 


130,922 


-204 


rook-8-8 


64 


40,320 


1,441,729 


-6,209 


match-9 


36 


945 


2,620 


-28 


match- 10 


45 


945 


9,496 


-1,216 


match- 11 


55 


10,395 


35,696 


-936 


match- 12 


66 


10,395 


140,152 


12,440 


match- 13 


78 


135,135 


568,503 


23,672 


mcgraph-7-2 


21 


217 


1,014,888 


-120 


nicgraph-8-2 


28 


504 


166,537,616 


-720 


mcgraph-9-2 


36 


1,143 


50,680,432,112 


-5,040 


randomv20gl00 


20 


100 


86,116 


-25 


randomv20g500 


20 


500 


227,792 


1,166 


randomv20gl000 


20 


1,000 


287,689 


-1,007 


randomv25gl00 


25 


100 


1,223,224 


-202 


randomv25g500 


25 


500 


3,628,979 


-3,815 


randomv25gl000 


25 


1,000 


5,368,430 


3,666 



Table 5. Characteristics of the examples used in Table [51 
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Frobby 


Frobby 








Example 


DBMS 


BCRT 


Sage 


Macaulay 2 


GAP 


rook-6-6 


0.01 


0.01 


1.04 


0.24 


0.13 


rook-7-7 


0.13 


0.14 


12.59 


3.37 


3.86 


rook-8-8 


2.43 


6.39 


>223.11* 


58.16 


>7200 


match- 9 


0.00 


0.00 


0.16 


0.13 


0.08 


match-10 


0.02 


0.01 


0.69 


0.29 


0.12 


match- 11 


0.21 


0.15 


2.90 


2.41 


6.97 


match- 12 


0.33 


0.47 


12,31 


8.32 


9.15 


m;Uch-l;] 


7.84 


11. 2() 


>7200 


101.49 


2401.58 


nicgraph-7-2 


0.00 


0.00 


322.94 


0.43 


22.41 


nicgraph-8-2 


0.03 


0.06 


>7200 


10.33 


>7200 


mcgraph-9-2 


0.40 


0.65 


>7200 


306.28 


>7200 


randomv20fl00 


0.00 


0.00 


11.37 


0.07 


0.42 


randomv20f500 


0.01 


0.01 


35.99 


0.43 


3.28 


randomv20fl000 


0.02 


0.02 


47.32 


0.72 


7.60 


randomv20fl00 


0.00 


0.00 


322.80 


0.37 


58.17 


randomv20f500 


0.02 


0.04 


>7200 


3.16 


592.70 


randomv20fl000 


0.02 


0.04 


>7200 


6.43 


>7200 



Sage reports taking 223.11s on rook-8-8, but the actual time was in 
excess of half an hour. The discrepancy persisted across several runs. 
We do not know how to explain the discrepancy since the times that 
Sage reports is usually in line with external measurements. 

Table 6. Comparison of Euler characteristic implementations. All times are in seconds. 



